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We study the structure and formation of the (j> mesic nuclei to investigate the in-medium 
modification of the <j!>-meson spectral function at finite density. We consider (p,(j>), (7,p) 
and (tt~ , n) reactions to produce a </>-meson inside the nucleus and evaluate the effects of its 
medium modifications to the reaction cross sections. We also estimate the consequences of 
the uncertainties of the K selfenergy in medium to the ^-nucleus interaction. We find that 
it may be possible to see a peak structure in the reaction spectra for the strong attractive 
potential cases. On the other hand, for strong absorptive interaction cases with relatively 
weak attractions, it is very difficult to observe clear peaks and we may need to know the 
spectrum shape in a wide energy region to deduce the properties of cj>. 

§1. Introduction 

Mesic atoms such as pionic- and kaonic-atoms are Coulomb assisted meson- 
nucleus bound systems and have been studied systematically for a long timeP The 
binding energies and widths of these bound states provide us unique and valuable 
information on the meson-nucleus interactions. Mesic nuclei, another kind of meson- 
nucleus systems, are meson-nucleus bound states mainly due to attractive strong 
interactions and have also been studied intensively in these daysP - ^ 1 In the con- 
temporary hadron- nuclear physics, mesic atoms and mesic nuclei are considered to 
be very interesting objects in the following two aspects. First, they are strongly in- 
teracting exotic many body systems which should be explored by nuclear physicists. 
Second, mesic atoms and mesic nuclei provide unique laboratories for the studies 
of hadron properties at finite density which are significantly important to explore 
various aspects of the symmetries of the strong interaction.®® 

In this letter, we report our studies on structure and formation of <f> mesic nuclei. 
From the study of the 4> meson properties in nucleus, we can intuitively expect to 
obtain new information on the OZI rule^ 1 and the ss component of the nucleon 
at finite density. The analysis of QCD sum rulesP and the data taken at KEKP 
suggested a 3 % mass reduction of (j) at normal nuclear density. On the other hand, 
in-medium properties of the eft meson have also been studied theoretically with a close 
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relation to K and K meson properties in medium because of the strong (ft KK 
coupling. The (ft meson selfenergy calculated in Refs. [T3"|). ll4p indicate a significantly 
smaller attractive potential for (ft. Since we can expect to obtain new information 
which is complementary to the invariant mass measurements,^ we investigate the 
(ft mesic nuclear states to study the (ft properties in medium. Especially, in the 
experiment}^ only a small fraction of the produced (ft meson decays inside the 
nucleus because of the long life time of the (ft meson. In the (ft mesic nucleus, all (ft 
mesons must stay inside the nucleus until their decay or absorption. Thus, we may 
have possibilities to get better information on (ft properties in nucleus. 

We show the calculated results of the bound states and formation spectra for 
some reactions in the following sections which will help to consider experimental 
feasibilities. 



§2. Optical potential and bound states 



We study the properties of the (ft mesic nuclei by solving the Klein-Gordon 
equation, 

[-V 2 + ^ + n4E,p(r))}(ft(f)=u 2 (ft(f) . (2-1) 

Here, p is the ^-nucleus reduced mass. We employ the empirical Woods-Saxon form 
for the nuclear density as, 

where R and a indicate radius and diffuseness parameters for the nucleus, and are 
fixed to be R = 1.1871 1 / 3 — 0.48 (fm) and a = 0.5 (fm), respectively. We solve the 
Klein-Gordon equation numerically following the method of Oset and SalcedoPS 
We show the (ft meson optical potential V^ pt (r,E) in Fig. [TJ which is obtained 

from the (ft meson selfenergy IJ^ in Ref. [T4"|) as V^ pt (r,E) = — II ^(E, p(r)) within 

the local density approximation. The selfenergy 77^, depends on the energy E of the 
(ft meson and on the density p. 

As shown in Fig. [TJ the theoretical potential based on 11^^ is a weak attractive 
potential, Re V^ pt (0,E = m^) ~ —7.5 MeV. The absorptive part of the potential 
Im V^p t (r, E) is also weak and has a similar strength as the real part. The theoretical 
model of Ref. 14) is based on the OZI rulep^ where the nuclear medium effects of 
the (ft meson are mainly from those of K and K through the strong (ftKK coupling. 
Hence, it seems natural that the ^-nucleus potential is weak. On the other hand, 
a 3% mass reduction is indicated in Refs. [TT]1 . [T2]) . To simulate it, we multiply 
by a factor the real part of the theoretical potential so that Re V^ pt (0,E = m^) 
scales to —30 MeV. We use both theoretical (shallow) and scaled (deep) potentials 
and compare the results to know the sensitivity of the observables to the potential 
depth. 

We calculate the binding energies B.E. and widths r of the bound states, which 
are related to the eigen energy u in Eq. (|2-1|) as u = p — B.E. — iT/2, by solving 
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Fig. 1. The <\> meson optical potentials as a function of the radial coordinate r for (f> - n B systems 
obtained from 77^ in Ref. I14p . The left and right panels show the real and imaginary parts. 
The solid, dashed, dotted, dotted-dashed lines indicate the potential strength for the <f> meson 
energies E - m = MeV, -10 MeV, -20 MeV and -30 MeV, respectively. 



the Klein-Gordon equation selfconsistently for the real part of the <f> meson energy 
to know the resonance energies in the reaction spectra observed in experiments. The 
calculated results are shown in Table H for four nuclei which are expected to be in 
the final states in the one proton pick-up reactions considered in the next section. 
We see that there exist several bound states of 4> in nuclei, however, the widths of 



Table I. Calculated binding energies and widths of <f> meson bound states in n B, 39 K, 123 In, and 207 T1. 
Widths do not include the 4> decay width in vacuum r^ ree = 4.26 MeV 23 in this Table. 



Shallow potential Deep potential 
B.E. [MeV] (r [MeV] ) B.E. [MeV] (r [MeV] ) 



n B 


none 


Is 


10.2(17.9) 


39 K 


none 


Is 


28.5(20.0) 






2p 


11.8(19.8) 


123 In 


Is 2.34(21.6) 


Is 


34.5(18.7) 






2s 


12.2(19.9) 






2p 


26.3(19.9) 


207rp| 


Is 3.73(22.5) 


Is 


40.1(17.5) 






2s 


24.4(20.2) 






2p 


27.1(20.3) 






3p 


13.1(20.4) 






3d 


27.5(19.8) 



the states are large and will make it difficult to observe any peak structure in the 
missing mass spectra. 
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§3. Formation spectra of (ft mesic nuclei 

We consider three types of proton pick-up reactions to form the <^>-mesic nucleus 
which are (7,p), (ir~ ,n) and (p, (ft). The elementary processes of these three reactions 
are 7 + p — > (ft + pp3 n~ + p — > (ft + np® and p + p — > (ft + 0™ respectively. We 
would like to make a few comments on the data in Ref. IT7|) . which have not been 
published yet. In the present calculation, we need to use the backward cross section 
of (ft production in the elementary process, which was evaluated here by a functional 
fit to the data shown in Ref. I17p . However there are large uncertainties in the value 
of the backward cross section that could affect the absolute value but not the shape 
of the spectrum. 

We first show in Fig. [2] the momentum transfer of the reactions as it is an 
important information of meson bound state formationPS'E] 1 \y e found that the 
momentum transfer of the (p, (ft) reaction is relatively smaller than the other two 
reactions. However, we also found that the momentum transfer is always larger than 
100 MeV/c in the energy range plotted in Fig. [2J 
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Fig. 2. Momentum transfer of the (-y,p) (solid curve), (ir~ ,n) (dashed curve) and (p,<f>) (dotted 
curve) reactions for the <f> mesic nucleus formation. Protons in the initial state are assumed to 
be in IP3/2 state in 12 C. 




We use the Green's function method to calculate the formation spectra of the 
(ft meson from nucleus^ 1 as in the cases of other mesic nucleus formationP^ \y e 
briefly explain the formalism for the (p, (ft) reaction as an example. The expected 

spectra of the (p, (ft) reactions ( — \ are evaluated from the nuclear response 



df2dE ( f > 

/ d(j \ ele 

function S(E) and the elementary cross section m the impulse approxima- 

tion as 

d 2 a \ /da\ ele 

1 — I * S(E) . (3-1) 



dQdE^J \dQJ 



The calculation of the nuclear response function with a complex potential is 
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Fig. 3. Calculated spectra for the 0-nucleus systems formation plotted as a function of the emitted 
4> meson energy in the 12 C(p, <f>) reaction at pp = 1.3 GeV/c. Total cross sections (upper panels) 
and conversion parts (lower panels) are calculated with the theoretical (shallow) optical potential 
(left panels) and the scaled (deep) optical potential (right panels). The dashed, dotted lines 
show the contributions from different proton hole states as indicated in the figure. The vertical 
dashed line indicates the <f> meson production threshold. 



formulated by Morimatsu and YazakP^ 1 as 

S(E) = --lmJ2 [ drdr'T ] f G{E-r,r')T f , (3-2) 

where the summation is taken over all possible final states. G(E; f, r') is the Green's 
function of (j) meson interacting in the nucleus and defined as, 

G(E;f,r') = (a\<f>(r) * ej>+(f')\a) , (3-3) 
E - H^ + te 

where a indicates the proton hole state and indicates the Hamiltonian of the 
(f> meson-nucleus system. The amplitude tj denotes the transition of the incident 
particle (p) to the nucleon-hole and the outgoing <j> meson, involving the nucleon-hole 
wavefunction ipj N and the distorted waves \i an d Xfi °f the projectile and ejectile. 
By taking the appropriate spin sum, the amplitude tj can be written as, 

r f (r) = X*f{m/2,m s K^) ® i'i^\jMXi{rl , (3-4) 

with the meson angular wavefunction Yuif) and the spin wavefunction £i/2 )ms of the 
ejectile. We assume the harmonic oscillator wavefunctions for ipj N with the empirical 
range parameter. 

The semi-exclusive spectra can be calculated by decomposing the response func- 
tion (|3-2p into the escape and conversion parts: S = S esc + S'con- This decomposition 
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can be done exactly by 




= E / drdr'r}(l + &V*(r, E))1mG (l + Vf pt (r,E)G)r f 



f J 



S C on(E) = --J2 I drdf'T}GhmV^ t (r,E)Gr f 



(3-5) 



where V3 t (r, E) is the (ft meson-nucleus optical potential given in the Hamiltonian. 
The conversion part is known to express the contributions of the (ft meson absorption 



We show the calculated cross sections of the 12 C(p, (ft) reaction proposed in 
Ref. [23]) for the formation of a eft - 11 B system in Fig. [3] for the theoretical (shallow) 
optical potential case in the left panels and for the scaled (deep) optical potential 
case in the right panels. The conversion part of the spectra are also shown in the 
lower panels in addition to the total spectra in the upper panels. The spectra shown 
in this article are folded spectra with a Gau/3ian distribution with (ft meson decay 
width in vacuum at E = m^. For the shallow optical potential case, since no bound 
state exists as shown in Table [U the spectrum has a smooth shape without any peak 
structures in the bound region. For the deep potential case, we find some enhance- 
ment of the spectra in the bound energy region. However, we do not observe any 
peak structures again for the deeper potential case even when there exists one bound 
state. We also find that the conversion contributions shown in lower panels, which 
correspond to the (p, (ft) spectra coincident with the particle emissions from (ft meson 
decay (absorption) in nucleus, resemble each other for both the shallow and deep 
potential cases, even though they have different constitutions of subcomponents. 
The conversion part of the spectra was defined in Eq. ()3-5p as a contribution of the 
absorptive interaction Im V^ t (r, E). One of the advantages of the (p, (ft) reaction is a 
possible background reduction by the conversion spectrum as described in Ref. [23]) . 
Thus, the conversion spectrum in Fig. [3] will be an important piece of information 
to observe the (ft meson properties in nucleus. 

To investigated the mass number dependence of the formation spectra, we show 
the calculated spectra in Fig. H]for 124 Sn target with the shallow and deep potential 
cases. We find that there appear some enhancements in the (ft meson bound energy 
region for the deep potential case. However, again we do not see any clear peaks 
corresponding to the bound states. As we have expected from the widths of the 
bound states, it is difficult to observe any peak structure in the standard missing 
mass spectra. 

In Figs. [5] and El we show the calculated results for different formation reactions, 
namely, 12 C(7,p) (Fig. [5]) and 12 C(7r~,n) (Fig. [6]). In both cases, the momentum 
transfers of the reactions are larger than for the p induced case and, hence, the 
spectra around threshold are suppressed in these reactions. 

We have checked the sensitivity of the optical potential vf pt (r,E) to the K~ 
selfenergy 77 K - in the nucleus. Since in the theoretical model of Ref. [H]), the medium 
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Fig. 4. Calculated spectra for the 0-nucleus systems formation plotted as a function of the emitted 
</> meson energy in the 124 Sn(p, <fi) reaction at pp = 1.3 GeV/c. We include the contributions 
from outer proton orbits, lgg/2 (dashed lines), 2p 1/ / 2 (dotted lines), and 2p 3 / 2 (dash-dotted lines). 
Total cross sections are calculated with the theoretical (shallow) optical potential (left panel) 
and the scaled (deep) optical potential (right panel). The vertical dashed line indiates the <j) 
meson production threshold. 



modification of the K + and K~ cloud around the 4> meson dominantly causes the 
modifications of the 4> properties, and since the K~ selfenergy Uk~ has a large 
uncertainty, it is interesting to change H K - and study the effects in 77^. We simply 
parametrize the modified K~ selfenergy 77^-- as 

77^- = a(Re 77^-) + 6(Im 77^-) , (3-6) 

where the 77^- in the right hand side is obtained by the chiral unitary model in 
Ref. rH|) . We consider the constant parameters a and b to scale the real and imaginary 
parts of n K ~. We show the calculated 77^ for four cases with a, b = 1 and/or 2 
in Fig. [71 We find that a certain ambiguity of 77^, due to that of 77^-- exists as 
we expected. However, the strength of the real part and imaginary part of TZ^ are 
correlated, namely, the strong attractive real part always appears with the strong 
absorptive imaginary part and so on. Thus, in any case, it will be difficult to see any 
clear structure well separated from the threshold and the quasi-elastic contribution 
in the formation cross section. 

Finally, we consider deep potentials to check the possibility to observe a peak 
structure. We scale again the theoretical potential in Ref. I14|) and consider deep 
real potential cases as Re 1^(0,75 = m^) = —40, — 70, —100 MeV keeping the 

Im Vf pt (r, E) to the theoretical value. The calculated results are shown in Fig. [HJ We 
find that we can observe a clear peak structure if the optical potential is as strongly 
attractive as |Re V* pt (0,E = m^)\ > 70 MeV. 

§4. Summary 

In this article we have studied the structure and formation of the (j) mesic nuclei 
to investigate the in-medium modifications of the (f) meson properties. We have 
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Fig. 5. Calculated spectra for the 0- nucleus systems formation plotted as a function of the emitted 
proton energy in the 12 C(7,p) reaction at p 7 = 2.7 GeV/c. The dashed, dotted lines show 
the contributions from s7, 2 and proton hole states, respectively. The vertical dashed line 
indicates the <f> meson production threshold. The theoretical (shallow) optical potential is used. 
The absolute value of this spectrum is still ambiguous. See detail in text. 



calculated the binding energies and the absorption widths of the <j) meson bound 
states. We have also studied the production of the eft meson in the nucleus and 
obtained numerical results for the formation spectra of the (p,(f>), (7,p) and (tt~,ti) 
reactions. Due to the large decay width of the 4> meson in nucleus, we found that 
the observation of a clear peak in the missing mass spectra will be difficult if the (f)- 
nucleus potential is weakly attractive as reported in Refs. fT3|) . [T^|) . A deep potential 
(Re V^p t (0, E = m^) ~ —30 MeV) corresponding to the 3 % mass reduction reported 
in Refs. [TT ]) .ll2 p is not enough to produce clear peaks in the reaction spectra for 
the absorptive potential strength evaluated in Ref. IT^|) . By scaling the strength 
of the real part of the potential, we find that we can expect clear peaks in the 
cross section due to the bound state formation for a strong attractive potential like 
| Re V^ pt {0,E = m^)\ > 70 MeV. We have also considered the sensitivity of the 
0-nucleus interaction to the uncertainties of the kaon selfenergies. 

As a conclusion, we find that it is difficult to observe clear signals of the </>■ 
mesic nuclei in the standard missing mass observation. However, the sub-threshold 
meson production in nucleus has the definite advantage to observe meson decay and 
absorption at finite density. Thus, we think that further theoretical investigations 
are very important to discover the better way to observe <p inside nucleus. 
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Fig. 6. Calculated spectra for the (^-nucleus systems formation plotted as a function of the emitted 
neutron energy in the 12 C(7r~,n) reaction at p n - = 2.0 GeV/c. The dashed, dotted lines show 
the contributions from sT^ and pT^ proton hole states, respectively. The vertical dashed line 
indicates the <j) meson production threshold. The theoretical (shallow) optical potential is used. 
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Fig. 7. Real part (left) and imaginary part (right) of the calculated ^-nucleus optical potential 
V^p t (r, E) = II^Etj,, p(r))/2/x are plotted as a function of the <j> meson energy at p(r) = po- 
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in text. 
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Fig. 8. Calculated spectra for the 0-nucleus systems formation plotted as function of the emitted 
0-meson energy in the 12 C(p, <f>) reaction at pp = 1.3 GeV/c. Total cross section (upper panels) 
and conversion parts (lower panels) are calculated with deep potentials obtained by scaling the 
theoretical potential as Re V£, t {0,E = m^) = -40, - 70, and -100 MeV while keeping 
the imaginary part unmodified. The dashed and dotted lines indicate the contributions from 
different sub components [IJ 1 ® as shown in the figure. 
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